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ABSTRACT 


The  coalescence  of  two  spheromaks  is  described  using 


the  quasi-ideal  model  of  reconnection.  Nonaxisymmetric 


effects  are  found  to  play  an  important  role  in  general. 
Possible  applications  of  repeated  coalescence  are  studied, 
including  heating  the  spheromak  plasma,  building  up  a  large 
spheromak  from  smalter  ones,  and  maintaining  a  spheromak  in 
steady  state.  The  analysis  suggests  that  coalescence 


experiments  could  also  provide  a  valuable  method  for 
studying  magnetic  field  relaxation. 
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I.  INTRODUCTION 


One  of  the  important  virtues  that  the  spheromak  shares  with  other 
compact  torus  devices  is  the  absence  of  a  toroidal  field  coil  linking 
the  plasma,  which  greatly  simplifies  the  engineering.  A  further 
consequence  of  this  geometry  is  the  possibility  of  moving  the  plasma 
along  the  external  vertical  field  (making  possible  a  moving  ring  field 
reversed  reactor*).  In  this  paper,  we  explore  yet  another  consequence 
of  this  geometry,  the  possibility  of  repeated  coalescence  of  separately 
formed  spheromaks.  Coalescence  of  two  partially  distinct  spheromaks  has 
in  fact  been  observed  on  the  University  of  Maryland  experiment, 

PS-1.2  Spheromak  coalescence  has  potential  applications  to  heating  the 
spheromak  plasma,  to  building  up  a  large  spheromak  from  smaller  ones 
(thus  making  it  easier  to  attain  fusion  conditions),  to  maintaining  a 
spheromak  In  steady-state,  and  to  the  experimental  study  of  magnetic 
field  relaxation. 

The  astron  experiment  at  Livermore  attempted  to  use  the  repeated 

3 

coalescence  of  e-layer  pulses  to  produce  a  field-reversed  e-layer. 
Coalescence  of  e-layer  pulses  was  observed,  although  field  reversal  was 
not  achieved  in  these  experiments.  Later  experiments  at  Cornell 
demonstrated  the  merging  of  two  electron  rings  below  field  reversal  to 
produce  a  field-reversed  electron  ring.^  Numerical  simulations  of  the 
astron  experiment^  also  looked  into  the  possibility  of  maintaining  a 
field  reversed  e-layer  by  repeated  injection  of  smaller  pulses.  These 
simulations,  assuming  axlsymmetry,  came  to  a  negative  conclusion. 

The  merging  of  ion  rings  has  been  proposed  as  a  method  of 
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maintaining  a  field  reversed  ion  ring,  and  axisymmetric  numerical 
simulations  are  in  progress. ^ 

The  recent  interest  in  spheromaks  has  sparked  several  proposals 
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concerning  spheromak  coalescence.  »  *  In  this  paper,  we  describe 
spheromak  coalescence  using  the  quasi-ideal  model  of 
reconnection. ^ ^  We  will  see  that  there  are  two  analytically 
soluble  limits:  one  in  whch  Kadomtsev's  theory  of  nonlinear  tearing 
modes can  be  applied,  and  one  in  which  Taylor's  relaxation  theory*  ■*» 
is  valid.  In  the  Kadomtsev  limit  there  is  an  infinite  set  of  conserved 
quantities  that  completely  determines  the  motion;  while  in  the  Taylor 
limit  only  one  of  these  quantities  remains  invariant,  the  system 
evolving  to  a  state  of  minimum  energy  subject  to  the  single  invariant. 

By  solving  in  both  limits  and  investigating  the  domains  of  validity  of 
both  approximations  we  develop  a  general  picture  of  spheromak 
coalscence.  Our  analysis  focuses  particularly  on  the  feasiblity  of 
maintaining  a  spheromak  in  steady  state  by  repeated  coalescence.  We 
also  calculate  the  heating  during  coalescence  for  some  particular  cases, 
and  look  at  the  feasibility  of  building  large  spheromaks  by  the  repeated 
coalescence  of  smaller  ones.  The  important  role  that  various  nonlinear, 
nonaxisymraetric,  tearing  processes  can  play  during  coalescence  will 
suggest  that  spheromak  coalescence  can  also  provide  a  valuable 
experimental  framework  for  studying  these  processes. 

An  axisymmetric  coalescence  is  shown  schematically  in  Fig.  1. 
Depending  on  the  force  exerted  on  the  spheromaks  by  the  externally 
applied  magnetic  field,  and  depending  on  the  force  due  to  the  spheromak 
fields  (the  spheromaks  attract  if  their  toroidal  currents  are  in  the 
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same  direction),  two  Initially  separated  spheromaks  can  approach  each 
(Jther  bv  moving  along  the  vertical  field.  As  the  spheromaks  merge,  the 
corresponding  flux  surfaces  reconnect.  After  the  process  has  gone  to 
completion,  the  two  spheromaks  have  combined  to  form  a  single  spheromak. 

Coalescence  has  been  observed  on  the  University  of  Maryland 
spheromak  experiment,  PS-1.4  Depending  on  the  configuration  of  the  bias 
magnetic  field,  the  initial  formation  stage  in  this  device  sometimes 
produces  two  separate  magnetic  islands,  which  subsequently  coalesce  on  a 
time  scale  short  compared  to  the  resistive  decay  time  of  the  plasma 
currents.  The  coalescence  involves  the  reconnection  of  a  substantial 
fraction  of  the  closed  flux  surfaces.  In  these  experiments,  the  initial 
formation  appears  to  he  well  described  by  an  axisymmetric,  resistive  MHD 
code. ^ 

Coalescence  of  tokamak  plasmas  is  also  possible  in  principle. 
Shafranov  has  proposed  heating  tokamak  plasmas  by  initially  forming  a 
number  of  separate  rings.  Inside  an  elongated  vacuum  chamber,  which 
subsequently  merge. His  analysis  assumes  that  the  reconnection 
process  is  axisymmetric,  and  that  the  plasma  is  incompressible.  The  use 
of  repeated  coalescence  to  maintain  the  toroidal  current  in  a  tokamak 
plasma  has  also  been  proposed;^  however,  the  fact  that  the  toroidal 
field  coll  must  link  both  the  main  and  refueling  plasma  rings  raises 
questions  concerning  the  practicality  of  repeated  coalescence  in  the 
case.  (Note  that  the  changing  poloidal  flux  through  the  hole  of  the 
refueling  ring  which  Initially  drives  its  toroidal  current  must  not  also 
thread  the  main  plasma  ring.)  The  analysis  of  Ref.  17  also  assumes 
axisymmetric  reconnection  and  plasma  incompressibility.  We  will  see 


that  we  cannot  rely  on  purely  axisymmetrlc  reconnection  to  maintain  the 
currents  in  a  spheromak. 

In  Sec.  II  we  look  at  purely  axisymmetrlc  reconnection,  and 
afterwards  include  nonaxlsymmetrlc  effects  In  Secs.  Ill  and  IV.  Sec. 

Ill  deals  with  two  step  coalescence  processes  consisting  of  an 
axisymmetrlc  reconnection  followed  by  a  nonaxlsymmetrlc  tearing  or 
double-tearing  mode.  We  apply  our  analysis  there  to  the  problem  of 
maintaining  spheromak  currents  in  steady-state,  and  calculate  the 
efficiency  of  the  refueling  process.  In  Sec.  IV  we  look  at  turbulent 
coalescence.  This  Is  found  to  be  of  particular  Interest  for  heating  the 
spheromak  plasma,  and  for  building  up  large  spheromaks  from  smaller 


ones. 
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II.  AXISYMMETRIC  RECONNECTION 


We  begin  by  looking  at  axisymmetric  reconnection,  and  afterwards 
include  nonaxisymmetric  effects  in  Secs.  Ill  and  IV.  The  assumption  of 
axisymmetry  is  a  reasonable  one  if  the  initial  and  final  spheromak 
equilibria  are  stable  to  nonaxisymmetric  modes.  We  will  see  that  this 
condition  is,  in  practice,  quite  restrictive. 

Having  assumed  axisymmetry,  our  flux  surface  correspond  to  level 
surfaces  of  the  poloidal  flux  function,  (r ) .  Flux  surfaces  which 
reconnect  must  have  the  same  value  of  ip.  We  use  the  quasi-ideal 
m<fdel^‘^*^  to  describe  the  reconnection  process.  The  plasma  is 
assumed  to  obey  ideal  MHD  except  in  a  narrow  region  about  the  x-point. 
When  flux  surfaces  reconnect,  the  enclosed  toroidal  fluxes  add.  Using 
the  fact  that  xty)  an  invariant  in  ideal  MHD,  we  get 

xf(¥>  -  xt(+)  +  x2<*)  »  0> 

where  x  denotes  the  toroidal  flux,  the  subscript  f  corresponds  to  the 
final  equilibrium  state,  and  the  subscripts  I  and  2  correspond  to  the 
initial  spheromaks.  Equation  (1)  applies,  of  course,  only  if  the  flux 
surface  ip  lies  Inside  both  spheromak  plasmas  initially.  We  take  the 
plasma,  hut  not  necessarily  the  plasma  current,  to  extend  to  the 
separatrix.  In  general,  the  initial  spheromaks  can  have  different 
values  of  ip  on  the  magnetic  axis,  ^  and  ip2*  If  lip^l  >  lip^l  > 
then  XfOlO  "Xi^)  f°r  l^l  ^  I'M  <  lipjl  •  The  value  of  }  on  the 
magnetic  axis  after  coalescence  is  ipj. 


Prom  Eq.  (1),  and  from  the  relation  q(i|i)  =  dx/<h|i  (where  q  Is  the 


safety  factor),  it  follows  that 

qf(t>  =  q^ifr)  +  •  ( 2 ) 

Either  Eq.  (1)  or  Eq.  (2)  can  be  viewed  as  specifying  an  Infinite  number 
of  constants  of  motion  for  axisymmetrlc  coalescence.  The  total  magnetic 
helicity,  K,  used  in  Taylor's  relaxation  theory,^*  ^  Is  one  of  these 
invariants, 


K=J  A  •  B  =  -2  /  x('l')  =*  2  /  i|iq(i|i)  di|i  ,  (3) 

v 

where  we  have  used  the  boundary  conditions  x  ■  0  on  the  magnetic  axis 
and  \|>  -  0  on  the  plasma  surface  to  Integrate  by  parts. 

The  poloidal  current  distribution  can  be  expressed  In  terms  of  q, 

I(i(i)  *  q(i|i)  di|i/dV  (<1/R2>^)"1  ,  (4) 

where  <  denotes  an  average  over  a  flux  surface,  V  Is  the  volume 
enclosed  by  the  flux  surface,  and  R  Is  the  radius  in  cylindrical 
coordinates.*®  Together  with  the  plasma  pressure,  I(i|i)  determines  the 
equilibrium  through  the  Grad-Shaf ranov  equation 


V%  -  R2p'(i|i)  +  II' (iji)  , 


7*  =  R  L1L+  3 

R  3R  R  3R  . 

3z 


where 
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Eliminating  I(i|i)  from  the  Grad-Shaf ranov  equation  in  terms  of  the  known 

we  obtain  a  QDE  which  determines  the  equilibrium.**  The  boundary 

condition  on  at  infinity  is  given  by  the  external  fields.  The 

plasma  f}(=  8itp/B  )  must  not  be  too  large  if  the  equilibrium  is  to  be 

1 9-2 1 

stable  to  Mercier  modes.  Since  the  poloidal  and  toroidal  betas  are 

comparable,  it  follows  that  the  spheromak  equilibrium  will  be  determined 

by  the  q  profile  alone  to  a  good  approximation. 

Now  we  can  look  more  closely  at  the  conditions  under  which  the 

assumption  of  axisymmetry  is  valid.  First  we  show  that  two  coalescing 

spheromaks  with  decreasing  q  profile  must  have  the  same  value  of  total 

poloidal  flux  to  preserve  axisymmetry.  In  Fig.  2  we  show  schematically 

what  happens  when  the  poloidal  fluxes  are  different.  The  final  q 

profile  is  not  monotonic,  and  is  therefore  unstable  to  ideal  modes  and 

to  double  tearing* J  modes.  Note  that  this  argument  is  specific  to 

decreasing  q  profiles.  For  initial  tokamak-like  q  profiles,  coalescence 

gives  a  q  profile  which  is  monotonic,  although  rapidly  varying  near  the 

i J /  corresponding  to  the  smaller  value  of  total  poloidal  flux. 

In  addition  to  requiring  equal  poloidal  fluxes,  we  must  also 

require  that  the  q's  on  the  magnetic  axis  for  the  two  coalescing 

spheromaks  sum  to  less  than  unity,  so  that  the  final  spheromak  will  be 
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stable  to  n  ■  1  internal  kink  modes.  If  their  peak  value  of  q  is  too 
small,  however,  the  initial  spheromaks  can  be  unstable  to 
nonaxisymraetric  ideal  modes.  Numerical  results  from  the  PEST  code 
Indicate  that  q  on  the  magnetic  axis  must  be  at  least  1/3.^ 

Now  we  turn  to  the  question  of  maintaining  a  spheromak  equilibrium 
in  steady-state  via  coalescence.  We  will  be  concerned  only  with 
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maintaining  the  q  profile  and  the  value  of  the  total  poloidal  flux. 
Throughout  this  paper  we  use  the  term  "refueling"  as  a  shorthand  for  the 
process  by  which  we  restore  the  spheromak  currents,  although  we  will  not 
be  concerned  with  restoring  the  plasma  density.  The  density  can  be 
increased  separately  through  the  injection  of  neutral  pellets.  Of 
course  we  must  require  that  coalescence  not  increase  the  plasma  density 
of  the  fusion  spheromak  above  its  initial  value.  (We  call  the  spheromak 
to  be  refueled  the  "fusion  spheromak".)  This  places  a  bound  on  the 
number  of  particles  allowed  in  the  refueling  spheromak, 

N2  S  Nl  0)(tM/tN)  * 

where  N2  an<*  Nl  are  the  numbers  of  particles  in  the  refueling  and  fusion 
spheromaks  respectively,  ip0  is  the  total  poloidal  flux  in  the  fusion 
spheromak,  A  ip  is  the  change  in  its  poloidal  flux  due  to  coalescence, 
and  tm  and  are,  respectively,  the  time  scales  for  flux  diffusion  and 
for  particle  loss.  As  we  have  mentioned,  the  plasma  $  must  not  be  too 
large,  so  q(^)  alone  will  determine  the  refueling  equilibrium  to  a  good 
approximation. 

Suppose  the  fusion  spheromak  initially  has  ip  *  on  its  magnetic 
axis.  After  some  time,  the  poloidal  flux  decays  to  a  value  ipQ  - 
To  restore  the  initial  value  of  poloidal  flux  by  purely  axisymmetric 
reconnection,  the  refueling  spheromak  must  have  ip  ■  on  its  magnetic 
axis;  but  we  have  seen  that  axisymoetry  is  not  preserved  unless  the 
poloidal  fluxes  of  the  coalescing  spheromaks  are  equal.  We  conclude 
that  it  is  impossible  to  maintain  a  spheromak  in  steady-state  by  purely 


axlsymmetrtc  coalescence.  That  conclusion  mav  appear  somewhat 
paradoxical,  since  we  can  use  Fq.  (2)  to  find  a  q  profile  for  the 
refueling  spheromak  which  restores  the  fusion  spheromak  to  its  Initial 
state. 

Again,  suppose  we  initially  have  a  spheromak  with  safety  factor  qQ 

on  the  magnetic  axis  and  total  poloidal  flux  iji  ,  as  shown  schematically 

o 

in  Fig.  3a.  After  some  small  time  has  elapsed,  both  these  quantities 
will  have  decayed  somewhat,  to  values  q  -  6q  and  —  6 \li  .  also  shown 
in  Fig.  3a.  (The  value  of  q  on  the  magnetic  axis  decreases  because  the 
current  profile  becomes  more  peaked  there.)  To  restore  the  initial  q 
profile  via  coalescence  requires  a  refueling  spheromak  whose  q(\|»)  is  the 
difference  of  the  two  q  profiles  of  Fig.  3a  (see  Fig.  3b).  The  required 
q  profile  is  discontinuous,  and  therefore  clearly  not  physically 
realizable.  The  obvious  thing  to  try  is  a  smoothed  version  of  this 
profile,  as  indicated  by  the  dashed  line  in  Fig.  3b.  Adding  this 
smoothed  profile  to  the  decayed  q  profile  of  the  original  spheromak,  we 
obtain  the  q(<|/)  shown  in  Fig.  3c.  The  equilibrium  corresponding  to  Fig. 

3c  is  unstable  to  nonaxi symmetric  modes  (e.g.  double  tearing  modes).  It 
is  clear  that  we  will  run  Into  the  same  problem  with  any  refueling 
spheromak  having  a  smooth  q(<p),  and  also  if  we  try  to  use  a  number  of 
different  refueling  spheromaks  in  succession. 

Although  the  equilibrium  corresponding  to  Fig.  3c  is  unstable,  we 
would  expect  the  resulting  turbulence  to  he  localized  near  ~  6^, 

and  to  result  in  a  smoothing  of  q(i|/).  This  suggests  a  possible  two  step 
process  for  restoring  spheromak  currents,  an  axisymmetric  coalescence  with 
a  refueling  spheromak  having  the  smoothed  q(<0  of  Fig.  3b,  followed  by  a 
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turbulent  smoothing  of  the  resulting  q  profile.  When  we  examine  the 

required  equilibrium  for  the  refueling  spheromak,  however,  this  does  not 

appear  to  he  a  very  attractive  scheme. 

We  solve  the  Grad-Shaf ranov  equation  analytically  in  the  limit 

that  6q  «  1  and  61I1  Ap  <<  1,  with  p  =  0,  I(ji)  given  by  Eq.  (4),  and 
o  o 

q(j,)  corresponding  to  the  smoothed  profile  in  Fig.  3b.  Our  approach  is 
to  match  solutions  which  are  valid  in  four  different  regions,  as 
follows: 

1)  Near  the  magnetic  axis  we  Taylor  expand  all  quantities 

in  . 

o 

2)  For  the  narrow  region  near  -  6\J>o  with  rapidly  varying  q  we 
derive  a  jump  condition. 

3)  This  is  the  region  where  q  is  small,  but  the  distance  from  the 
magnetic  axis  is  still  much  less  than  the  major  radius  of  the  spheromak. 

4)  Far  from  the  magnetic  axis  we  can  use  the  elliptic  function 

solution  for  the  fields  due  to  a  circular  loop  of  current. 

Further  details  are  given  in  the  Appendix.  We  find  that  the  equilibrium 

has  a  hollow  current  profile  and  a  large  flux  hole  (that  is,  a  large 

current  free  region  near  the  separatrix).  The  low  value  of  q  on  the 

magnetic  axis  (q  =  5q  there)  may  make  stability  hard  to  achieve. 

o 

Using  nonaxisymmetric  effects  to  produce  poloidal  flux  from 
toroidal  flux,  we  can  find  considerably  more  attractive  refueling 


schemes 
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III.  NONAXISYMMETRIC  EFFECTS  DUE  TO  SLOWLY  GROWING  MODES 


If,  having  assumed  that  the  coalescence  process  is  purely 

axisymmetric,  we  find  that  the  resulting  spheromak  is  unstable  to 

nonaxisyrametric  modes,  we  must  conclude  that  our  assumption  of 

axisymmetry  has  broken  down  at  some  point.  When  the  unstable  modes  grow 

sufficiently  slowly,  however,  we  can  still  treat  the  coalescence  as  a 

two  stage  process:  an  axisymmetric  reconnection  followed  by  a 

nonaxisymmetric  tearing  mode.  This  treatment  is  valid  if  the  e-folding 

time  of  the  fastest  growing  mode  is  long  compared  to  the  time  scale  on 

which  the  axisymmetric  reconnection  occurs.  For  example,  for  the  n  =  1, 

m  *  1  mode,  which  will  be  of  particular  interest  to  us  in  this  section, 

the  growth  rate  is  small  if  the  maximum  value  of  q  is  not  much  larger 

than  one.  For  an  equilibrium  having  q  *  1.3  on  the  magnetic  axis, 

Gautier  et  al.  found  a  growth  rate  of  .03  where  is  the  Alfven 

time.  By  comparison,  in  the  experiment  reported  in  Ref.  2,  two 

Islands  containing  about  25%  of  the  total  poloidal  flux  coalesced  in  a 

time  not  much  longer  than  the  Alfven  time.  In  a  numerical  study  of  the 

coalescence  of  two-dimensional  islands,  Pritchett  and  Wu  found  that  the 

reconnection  time  for  this  type  of  forced  reconnection  is  relatively 

2  5 

insensitive  to  the  value  of  the  resistivity. 

The  tearing  modes  we  will  be  concerned  with  in  this  section  will  be 
localized  sufficiently  close  to  the  magnetic  axis  that  helical  symmetry 
can  be  assumed.  The  flux  surfaces  are  determined  by  the  helical  flux 
function,  t  =  q^  -  y  ,  where  q8  is  the  value  of  q  at  the  separatrix. 
Flux  surfaces  which  reconnect  must  have  the  same  value  of  ♦ . 


We  describe  the  nonlinear  evolution  of  the  tearing  mode  using 


Kadomtsev's  theory.^0  Again,  the  plasma  Is  assumed  to  obey  Ideal  MHD 
except  In  a  narrow  region  about  the  separatrix.  Those  flux  surfaces 
which  initially  have  the  same  value  of  *  (hut  different  values  of  x)  all 
reconnect,  so  that  at  the  end,  ♦  uniquely  labels  the  flux  surfaces. 

When  flux  surfaces  reconnect  the  enclosed  toroidal  fluxes  must  be 
subtracted, 


Xf<*>  “  *  X2^)  * 

where  the  subscripts  1  and  7  distinguish  the  two  values  of  the  initially 
double-valued  function  an<1  is  the  single  valued  function 

which  results  after  the  reconnection  process  has  gone  to  completion. 

Kadomtsev's  theory  has  an  infinite  set  of  conserved  quantities.  It 

was  recognized  by  Bhattachar jee  et  al.  that  the  invariants  of 

f  3 

Kadomtsev's  theory  could  be  expressed  in  the  form  J  f(*)  A  •  B  d  x  , 

V 

where  V  denotes  an  integral  over  the  entire  volume  of  the  plasma  and  f 
is  an  arbitrary  function,  depending  on  the  helicity  of  the  tearing  mode 
we  get  different  helical  flux  functions,  ♦,  and  thus  different  sets  of 
conserved  quantities.  If  we  set  f ($ )  =  1,  we  get  the  total  magnetic 
helicity  employed  by  Taylor.  This  is  the  only  one  of  the  invariants 
which  is  preserved  by  a  mode  of  any  helicity.  Recall  that  this  quantity 
is  also  conserved  during  the  axisymmetric  reconnection. 

We  look  first  at  the  case  where  i|i  on  the  magnetic  axis  is  the  same 
for  both  coalescing  spheromaks,  and  where  the  q 's  on  the  axis  sum  to  a 
value  slightly  above  one.  This  will  be  of  particular  interest  for  the 
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refueling  problem.  Figure  4  shows  the  q  profiles  before  and  after  the 
axisymmetric  stage  of  coalescence.  Because  the  unstable  mode  Is 
localized  near  the  magnetic  axis,  we  can  treat  It  by  Taylor  expanding 
the  equilibrium  quantities  about  the  axis.  (We  also  make  use  here  of 
the  fact  that  the  Grad-Shaf ranov  equation  is  nonsingular  at  the  magnetic 
axis.)  For  the  equilibrium  after  the  first  stage  of  the  coalescence,  we 
write 


$  “  4> 

o 


+  ax  +  bx 


2 

t 


(6) 


where  $  Is  the  value  of  *  on  the  magnetic  axis,  and  a  and  b  are 
o 

constants  which  determine  the  equilibrium  in  the  neighborhood  of  the 
axis.  Expressing  ip  In  terms  of  *  and  x*  we  8et 

,),  «  *o  +  (a  +  1)  X  +  bx2  ,  (7) 

so  that  ♦  ■  ^  -  6ipQ.  Taking  the  conventions  x  >  0  and  q  ■  dx/<kp  >  0, 
we  conclude  from  ip  -  0  at  the  separatrlx  that  <  0.  From  Eq.  (7),  we 
calculate 


q  ■  1/(1  +  a  +  2bx)  «  (8) 

Comparing  this  result  with  the  q  profile  of  Fig.  4b  we  see  that  a  <  0 
and  b  >  0. 


Figure  5  shows  what  #(x)  looks  like  before  the  tearing  mode  begins 
to  grow  (solid  line).  The  minimum  value  of  •  lies  on  the  n  ■  1,  m  ■  1 


rational  surface,  which  corresponds  to  x  =  X  =  “  a/(2b)  .  The  value  of 

8 

2 

$  there  is4g=4>o-a  /(4b)  .  In  terms  of  4g  and  xg  we  can  rewrite 

2 

Eq.  (6)  as  4  -  4<,  =  b(x  -  Xg)  »  or 

X  =  xs  ±  ( (*  -  *s)/bl^2  .  (9) 

The  tearing  mode  leads  to  reconnection  of  flux  surfaces  having  the 

same  value  of  4.  From  Eqs.  (9)  and  (5)  we  see  that  after  reconnection 

1/2 

we  have  x  =  2  [(4  -  4-)/bl  ,  or 

4  =  b  x2/4  +  4S  ,  (10) 

which  is  shown  as  the  dashed  line  in  Fig.  5.  Only  those  flux  surfaces 
having  x  <  2  x$  are  affected  by  this  second  stage  of  reconnection.  The 
final  q  profile  in  this  region  is 

q  -  1/(1  +  b  x/2)  .  (11) 

Note  that  at  the  boundary  of  this  region,  Eq.  (11)  gives 
q  +  1/(1  -  a/2)  while  Fq.  (8)  gives  q  +  1/(1  -  a)  .  We  anticipate  that 
this  discontinuity,  which  is  characteristic  of  Kadomtsev's  theory,  is 
smoothed  hy  further  turbulence.  On  the  magnetic  axis  q  is  equal  to  one. 

Now  we  apply  these  results  to  the  refueling  problem.  To  do  so,  we 
require  that  the  two  stage  reconnection  process  we  have  just  described 
restore  the  total  poloidal  flux  to  ip0,  and  that  it  approximately  restore 
the  q  profile  to  that  of  the  fusion  spheromak  at  t  ■  0.  Letting  qQ0|>) 


r 


be  the  q  of  the  fusion  spheromak  at  t  -  0,  which  Is  assumed  to 

have  q  0|*  )  *•  1  to  maximize  the  allowable  beta,  our  condition  on  q  It 
o  o 

1  -  3a/4  -  [q  I  l-1  .  M2' 

°  X  "  2Xs 

This  gives 


b  -  -  (4/3)  q'  (\|>  )  , 

o  o 


where 


qo(,*'o)  =  ijrqo(*)l^o  • 

The  condition  on  the  poloidal  flux  then  gives 

a2  -  (16/3)  q;(*0>«4>0  . 

The  q(<i)  for  the  refueling  spheromak  is 

qj  *  Sq  +  ^2  +  (5/3)q^(tf»o)(t|i-^o)  (13a) 


for 


<  (4/5)|36tJ>  /q'(i  )| 
o  o  o  o 


(13b) 


and 


J 


"  Sq(ip) 


(14) 


for  larger  values  of 

We  can  solve  analytically  for  the  refueling  equilibrium  specified 

by  Eqs.  (13)  and  (14)  In  the  limit  of  small  Sq  and  6i|>  by  the  same 

o 

methods  used  to  solve  the  equilibrium  of  the  previous  section.  This  Is 
discussed  In  the  appendix.  The  equilibrium  is  now  considerably  more 
Interesting,  with  no  hollow  current  profile,  a  smaller  flux  hole,  and  a 
larger  value  of  q  on  axis.  The  condition  that  q 2 Cip Q)  >  .3,  required 
for  stability, ^  Is  satisfied  for  6t|i  /i|>  >  .02.  In  obtaining  this 

estimate  we  have  approximated  *  q  (i|>  )/il>Q. 

We  calculate  the  efficiency  of  our  process  for  maintaining  the 
spheromak  currents  in  the  limit  that  the  fusion  and  refueling  spheromaks 
both  have  a  large  flux  hole.  That  is,  we  determine  the  total  magnetic 
energy  lost  during  coalescence  as  a  fraction  of  the  energy  of  the 
refueling  spheromak.  This  lost  magnetic  energy  Is  dissipated  by  ohmic 
heating  near  the  x-polnt. 

The  magnetic  energy  of  each  of  the  (large  aspect  ratio)  spheromaks 
can  be  written  as 


W  -  i  LI2  +  Mil  ,  (15) 

L  e 

where  L  Is  the  self  inductance  of  the  spheromak  currents,  M  is  the 
mutual  inductance  between  the  spheromak  currents  and  the  external  field 
coil,  I  Is  the  total  toroidal  current  of  the  spheromak,  and  Ig  is  the 
total  current  through  the  external  field  coll.  The  inductances  are 
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related  to  the  flux  through  the  spheroraak  hole  by  the  relations 


MI0  *  ^  /c, 

e  e 


(16) 


and 


LI  ■  \pg/c  (17) 

where  U»  Is  the  external  flux  and  ib  Is  the  self  flux.  The  previously 
Te  s 

defined  pololdal  flux  function  Is  proportional  to  the  total  flux, 

e  s’ 

where  R  is  the  malor  radius.  We  define  ¥  =  2irRi|i  .  The  self  inductance 

o  o 

corresponding  to  the  refueling  equilibrium  is 

L£  "  (4xR/c^)Un(8/5qm)  -  2)  , 

where  6q  Is  the  value  of  6q  on  the  magnetic  axis.  To  estimate  the  self 
tn 

Inductance  of  the  fusion  spheromak  we  use  the  standard  expression  for 
the  inductance  of  a  large  aspect  ratio  circular  torus, 

Lt  -  (4vR/c2)Un(8R/a)  -  2  +  1^2], 

where  a  Is  the  minor  radius,  and  >  1/2  is  the  internal  Inductance  per 


unit  length 
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Equilibrium  force  balance  for  a  large  aspect  ratio  spheromak  gives 
2irRBvI/c  *>  (!/2)I2dL/dR  -  I2L/(2R)  ,  (18) 

where  Bv  Is  the  externally  Imposed  vertical  field.  From  Eq.  (18)  we  get 

2 

¥  =*  -(l/4)¥  and  Mil  *  -(1/4)LI  .  It  follows  that  the  magnetic  energy 

e  s  e 

of  each  spheromak  can  be  expressed  as 

W  -  (4/a)  ¥2/(c2L)  .  (19) 

o 

Using  Eq.  (19)  In  conjunction  with  Eq.  (18),  we  find  that  when  the  total 
poloidal  flux  of  the  fusion  spheromak  decays  by  an  amount  6¥q,  the 
magnetic  energy  of  this  spheromak  decays  by  an  amount 

5W.  =  (2/3)  ¥  6¥  /(L  c2)  .  (20) 

l  o  o  l 

If  our  refueling  process  were  perfectly  efficient,  the  magnetic  energy 
of  the  refueling  spheromak  would  .just  be  given  by  Eq.  (20).  For  the 
refueling  process  we  have  described,  the  efficiency  Is  given  by 

SWj/V^  -  3(6iPQ/to)(L2/L1)  -  3(5<PoA|/o)Un6qm/*n(a1/R1)]  .  (21) 

For  example,  for  5<|/  A| i  »  $q  a  .1  , 


A  substantial  fraction  of  the  magnetic  energy  of  the  refueling 
spheromak  gets  transformed  into  thermal  energy  during  coalescence. 

There  are  two  other  possible  sources  of  thermal  energy,  the  kinetic 
energy  of  the  spheromaks  due  to  their  relative  velocity  on  impact,  and 
adiabatic  heating.  There  could  also,  in  principle,  be  adiabatic 
cooling.  We  will  have  more  to  say  about  these  effects  in  the  next 
section,  whe'v  we  will  look  more  closely  into  the  possibility  of  using 
coalescence  to  heat  the  spheromak  plasma. 

inus  f?»r  in  this  section  we  have  assumed  that  the  coalescing 

sober viuj.'.  ?.  have  an  equal  amount  of  poloidal  flux.  Since  the  magnetic 

2 

energy  of  the  refueling  spheromak  is  proportional  to  it  clearly 

would  be  advantageous  if  the  poloidal  flux  of  this  spheromak  could  be 
smaller.  If  the  coalescing  spheromaks  have  different  amounts  of 
poloidal  flux,  axisymmetric  reconnection  gives  a  q  profile  with  a  local 
maximum  (as  in  Fig.  6).  The  corresponding  equilibrium  is  unstable  to  a 
double-tearing  mode.  Assuming  that  the  growth  rate  of  the  double¬ 
tearing  mode  is  sufficiently  small,  our  two  step  coalescence  process  now 
consists  of  an  axisymmetric  reconnection  followed  by  a  nonaxisymmetric 
double-tearing  mode.  We  use  Kadomtsev's  theory  of  the  nonlinear 
evolution  of  the  double-tearing  mode. 

The  n  *  1,  m  *  1  helical  flux  function  corresponding  to  the  q 
profile  of  Fig.  6  is  shown  schematically  in  Fig.  7  (dashed  line).  There 
are  two  q  *  1  rational  surfaces,  correponding  to  the  two  local  extrema 
in  4  (a  local  minimum  and  a  local  maximum).  We  again  Taylor  expand  4  as 
a  function  of  x»  expanding  about  x  “  0  for  X  ^  Xsi*  an<*  ab°ut  X  "  Xsi 
for  x  >  Xsi*  We  exPre88  the  coefficient  directly  in  terms  of  the  q 


profile  at  this  Intermediate  stage.  For  x  <  XS1»  we 


6q(4>  -  6t()  ) 

o  q,  (ty  -  64)  )  * 
1  To  o 


1  *\<*o  ~  6*Q>  ^2 

2  f<>l(*o  '  6* o}1 3 


with 


4  a  4>  —  64) 

0  0  o 


and 


Xsi  *  (*2  "  Vql(*o  '  6*0}  +  (1/2)(*2  "  *o)2q'l(*o  ‘  6V 


For  x  >  Xsi  ic  *s  convenient  to  define 


qm=  +  ^o  ’  V 


and 


dq. 


d4» 


4»  -  64) 

To  To 


dq , 

_ C 

d4* 


4)  -  d4» 

ro  o 


We  then  find,  for  x  >  Xsi» 


4*4.+  - 

1  <1. 


m  <x  -  x<5i> - ?  (x  -  xsl)2  » 


m 


‘S1  2q3 

m 


with  4j  "  ^2  ~  Xri  *  l°cal  minimum  of  Eq.  (25)  corresponds 


(22a) 


(22b) 


.(23) 


(24a) 


(24b) 


(25) 

to  the 


second  q  ■  1  rational  surface,  at 


(26) 


XS2  ~  XS1 


+  qO  -  q_)/q' 

m  Itl  m 


» 


4,  =  4,  +  (l/2)q  (1  -  q  )2/q' 
£.  i  m  mm 


(27) 


In  Kadomtsev's  theory,  all  flux  surfaces  having  the  same  value  of  4 

reconnect.  The  value  of  4  on  the  magnetic  axis  is  affected  by  this 

process  only  if  |4  _  |  >  |4  |  *  fib  -  64  |.  When  this  condition  is 
z  O  o  o 

satisfied,  the  final  value  of  4  on  the  axis  is  just  *2*  To  restore  the 
poloidal  flux  of  the  fusion  spheromak,  we  must  have  4j  =  4^,  or 


qraU  -  q  )2/q' 

Tn  m  m 


(4*2  “  -  «40)  +  26<Po  -  2(ip ^  ~  *Q)  «q 


(28) 


From  Eqs.  (22)  -  (27),  we  determine  what  x(®^  must  be  after  all 
corresponding  flux  surfaces  have  reconnected, 

1  1 1I 

!2[2(4  -  4,)q  /q  '  ]  for  4,  <  $  <  $ 

(29) 

2[2(4  -  *2)qm/qm]1/2  +  '  *0*>i<*0  "  for  *\  >  ♦  >  *Q  • 

The  flux  surfaces  with  4  >  4j  are  unaffected  by  the  n  *  1,  m  »  1  double¬ 
tearing  mode.  The  final  n  ■  1,  m  ■  1  helical  flux  function  is  shown 
schematically  in  Fig.  (7)  as  a  dashed  line.  The  corresponding  q(ip)  is 

not  monotonic  in  the  neighborhood  of  4  ■  4  .  We  could  expect  that 

o 

further  turbulence  will  smooth  out  this  q  profile. 

Equation  (28)  shows  that  as  4 ^  ~  •  increases,  so  does  the  required 
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value  of  qffl  -  1.  When  qm  -  1  Is  not  small,  the  growth  rate  of  the  n  * 

1,  m  ■  1  nonaxisymmetric  mode  Is  not  small,  and  we  can  no  longer  expect 
the  coalescence  to  proceed  In  two  separate  stages.  That  is,  we  can  no 
longer  make  any  symmetry  assumptions  about  the  coalescence  process. 
Similarly,  our  assumption  of  helical  symmetry  for  the  nonlinear  tearing 
modes  breaks  down  if  reconnection  occurs  sufficiently  far  from  the 
magnetic  axis  that  toroidal  effects  are  important,  or  if  tearing  modes 
of  different  helicity  grow  to  large  amplitude  simultaneously  and  begin 
to  overlap.  In  these  cases,  we  introduce  an  increased  amount  of 
turbulence  into  the  fusion  spheromak,  something  that  is  clearly 
undesirable  in  any  refueling  scheme.  Turbulent  coalescence  is,  however, 
of  interest  for  heating,  and  for  building  up  larger  spheromaks  from 
smaller  ones.  We  describe  turbulent  coalescence  in  the  next  section. 


A 
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IV.  TURBULENT  COALESCENCE 


Thus  far  we  have  been  assuming  that  at  each  stage  we  have  either 
axial  symmetry  or  helical  symmetry.  Within  the  quasi-ideal  model,  there 
is  an  infinite  set  of  constants  of  motion  associated  with  each  of  these 
symmetries.  It  is  these  invariants,  in  essence,  that  we  have  used  to 
determine  the  final  equilibrium  state  resulting  from  the  coalescence 
process.  Recall,  too,  that  the  total  magnetic  helicity  is  the  only 
quantity  which  is  conserved  by  a  tearing  mode  of  arbitrary  helical 
symmetry,  and  that  it  is  also  an  invariant  under  axial  symmetry. 

Our  symmetry  assumptions  can  break  down  for  a  variety  of  reasons, 

as  discussed  at  the  end  of  the  previous  section.  When  the  assumption  of 

symmetry  is  not  valid,  we  expect  that  there  will  be  fewer  invariants, 

and  that  the  plasma  will  assume  a  state  of  minimum  energy  subject  to 

whatever  set  of  invariants  is  appropriate.  As  long  as  the  quasi-ideal 

model  remains  valid,  the  total  magnetic  helicity  will  be  conserved.  In 

the  region  where  ideal  MHD  is  obeyed,  away  from  the  separatrix,  the 

local  magnetic  helicity  is  conserved.  Regardless  of  the  shape  of  the 

separatrix,  reconnection  of  two  flux  surfaces  clearly  does  not  affect 

the  total  magnetic  helicity.  In  the  limit  that  the  reconnection  process 

becomes  extremely  turbulent,  we  expect  that  only  the  total  magnetic 

helicity  is  conserved.  The  final  state  is  that  having  minimum  energy 

subject  to  the  conservation  of  magnetic  helicity.  This  is  the 

hypothesis  that  Taylor  has  used  successfully  to  describe  the  initial 

1 

formation  of  the  reversed  field  pinch.  We  use  Taylor's  theory  in  this 
section  to  describe  turbulent  coalescence. 
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For  our  coalescing  spheromaks,  we  specialize  to  the  spherical 
minimum  energy  equilibria  of  Rosenbluth  and  Bussac,  which  have  a 
uniform  external  magnetic  field  far  from  the  spheromak.  In  spherical 
coordinates,  the  poloidal  flux  function  in  the  plasma  is 


4>(r,6) 


sin^Q  jj(Xr)  /X 


(29) 


with 


V  x  B  =  XB  ,  (30) 

where  X  is  a  constant.  Although  this  equilibrium  is  unstable  to  a 
variety  of  modes  in  the  absence  of  a  conducting  wall,26  the  modes  can  be 
stabilized  by  a  combination  of  modifications  to  the  equilibrium  and  a 
relatively  distant  external  conducting  wall,2**^*2®*2^  neither  of  which 
should  strongly  affect  our  conclusions  regarding  turbulent  coalescence. 

The  equilibrium  solution,  Eq.  (29),  has  a  free  parameter,  X,  which 
is  a  function  of  the  total  magnetic  helicity,  K.  The  radius  of  the 
spheromak  is  proportional  to  while  the  total  poloidal  flux  is 

proportional  to  K^2,  and  the  total  magnetic  energy 
■  (l/2c)  /  J  •  A  d^x  *  XK/8x  is  proportional  to 

Now  suppose  that  we  coalesce  two  spheromaks  of  equal  size,  each 
having  a  total  magnetic  helicity  of  KQ.  The  final  equilibrium  must  have 
total  magnetic  helicity  2KQ.  If  the  turbulence  level  during  coalescence 
is  sufficiently  large  that  Taylor's  hypothesis  is  valid,  the  final 
equilibrium  after  coalescence  will  again  satisfy  Eq.  (30)  in  the  plasma. 


with  \  determined  by  K  =  2KQ.  With  plasma  out  to  the  separatrix,  the 
final  equilibrium  in  the  uniform  external  field  again  belongs  to  the 
class  of  spherical  Rosenbluth-Bussac  equilibria,  with  i(t(r,0)  in  the 
plasma  of  the  form  of  Eq.  (29).  The  total  poloidal  flux  in  the  final 
spheromak  is  40%  larger  than  in  each  of  the  initial  spheromaks,  while 
the  total  magnetic  energy  is  1.7  times  what  it  was  in  each  initial 
spheromak.  Fifteen  percent  of  the  total  magnetic  energy  has  been 
ohmically  dissipated. 

The  substantial  increase  in  poloidal  flux  due  to  Taylor  relaxation 
suggests  that,  in  coalescence  experiments,  measurement  of  the  total 
poloidal  flux  by  itself  would  give  a  useful  indication  of  what  is 
happening.  Recall  that  axisymmetric  coalescence  preserves  the  maximum 
value  of  \Ji.  (Due  to  the  presence  of  finite  resistivity,  i|>  actually 
decays  slightly.)  Nonaxisymmetrlc  tearing  modes  prod".;;*  poloidal  flux, 
increasing  somewhat  the  maximum  value  of  \|>. 

Next,  let  one  of  the  coalescing  spheromaks  be  very  much  smaller 
than  the  other.  Let  the  total  magnetic  helicities  be  Kj  and  K2,  with  K2 
*  eKj.  The  fraction  of  the  magnetic  energy  of  the  second  spheromak  that 
gets  dissipated  by  ohmic  heating  is  then 

AW/W2  -  1  -  (3/4)  e1/4  . 

In  terms  of  the  initial  radii  of  the  two  spheromaks,  this  can  be 
rewritten  as 


A W/W.  -  1  -  (3/4)  R./R, 
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If  Che  second  spheromak  is  much  smaller  than  the  first,  almost  all  of 
its  magnetic  energy  gets  transformed  into  thermal  energy  during 
coalescence*  It  appears  that  this  could  be  an  effective  means  of 
heating  the  spheromak  plasma.  This  effect  also  places  a  limit  on  the 
building  up  of  larger  spheromaks  by  repeated  coalescence  of  smaller 
ones. 

Suppose,  by  repeated  coalescence,  we  merge  n  spheromaks  of  equal 
size.  The  total  magnetic  energy  of  the  final  spheromak  is  n-*^  times 
the  sum  of  the  energies  of  the  initial  spheromaks.  The  radius  of  the 
spheromak  is  increased  by  a  factor  n^^. 

To  determine  the  change  in  thermal  energy  of  the  spheromak  plasma 

during  coalescence  we  should  include  also  the  contribution  due  to  the 

kinetic  energy  the  spheromaks  on  impact,  and  that  due  to  adiabatic 

heating  (or  cooling).  We  examine  the  kinetic  energy  contribution 

first.  Suppose  first  that  we  initially  have  an  equilibrium  state  with 

two  separated  spheromaks.  The  external  field  is  then  modified  in  some 

way  to  allow  the  spheromaks  to  coalesce.  We  show,  following 

Shaf ranov, ^  that  the  kinetic  energy  in  this  case  is  small  unless  the 

initial  separation  between  the  spheromaks  is  quite  large.  Let  At  be  the 

time  it  takes  to  make  the  necessary  changes  in  the  external  field,  and 

let  L  be  the  initial  separation  distance  of  the  two  spheromaks.  The 

kinetic  energy  density  on  impact  is  of  the  order 
7  2  2  2  2 

pL'/At  ■  (t./At)  (L/2R)  (B  /4x ) ,  where  p  is  the  local  mass  density  of 
A 

the  plasma,  t^  Is  the  Alfven  time,  R  is  the  major  radius  of  the 
spheromak,  and  B  is  the  local  magnetic  field.  In  general  t^  «  At,  and 
the  kinetic  energy  is  a  very  small  fraction  of  the  magnetic  field  energy 
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unless  L  »  R.  The  Initial  state  need  not  be  an  equilibrium  (In  fact, 
in  the  PS-1  experiment,  the  initial  state  with  separated  islands  is  not 
an  equilibrium);  nevertheless,  the  contribution  of  the  kinetic  energy  to 
the  energy  balance  is  not  significant  unless  the  velocity  of  the 
spheromaks  is  a  substantial  fraction  of  the  Alfven  velocity. 

We  turn  finally  to  the  question  of  adiabatic  heating  and  cooling. 
For  the  spherical  Rosenbluth-Bussac  equilibria,  r^  scales  like 
After  coalescence  the  plasma  occupies  a  slightly  smaller  total  volume 
than  it  does  initially,  so  there  is  a  slight  heating  effect.  For 
coalescence  of  two  equal  spheromaks,  for  example,  the  total  volume 
occupied  goes  down  by  about  15%.  The  effect  of  adiabatic  heating  is 
clearly  small  compared  to  that  of  ohmic  heating  near  the  x-point. 


V.  CONCLUSIONS 


Using  the  quasi-ideal  model  of  magnetic  reconnection,  we  have  been 
able  to  describe  spheromak  coalescence  in  two  analytically  soluble 
limits:  one  in  which  Kadomtsev's  theory  applies,  and  one  in  which 
Taylor's  theory  applies.  This  has  enabled  us  to  develop  a  general 
picture  of  spheromak  coalescence.  We  have  found  that  a  substantial 
fraction  of  the  magnetic  energy  is  transformed  into  thermal  energy 
during  coalescence  due  to  ohmic  heating  near  the  x-point.  Repeated 
coalescence  with  smaller  spheromaks  could  provide  a  particularly 
effective  means  of  heating  a  spheromak  plasma.  Repeated  coalescence 
could  also  he  used  to  build  a  large  spheromak  from  smaller  ones;  but  the 
relatively  large  amount  of  magnetic  energy  dissipated  when  a  small 
spheromak  coalesces  with  a  much  larger  one  places  limits  on  the  size  of 
a  spheromak  that  can  be  practically  formed  in  this  way.  We  have  also 
found  that  spheromak  currents  can  be  maintained  via  repeated  coalescence 
with  "refueling  spheromaks"  having  the  appropriate  q  profile  and  total 
poloidal  flux.  Nonaxisymmetric  effects  play  a  key  role  in  this 
process.  Finally,  our  analysis  suggests  that  coalescence  experiments 
could  lead  to  a  beter  understanding  of  the  Taylor  and  Kadomtsev  models 


in  nonlinear  MHD 
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APPENOIX 


In  this  appendix  we  obtain  analytic  solutions  for  our  refueling 
equilibria  in  the  limit  of  small  5q  and  6i|u  We  look  first  at  the 
equilibrium  discussed  at  the  end  of  section  II,  following  the  approach 
outlined  there  to  obtain  our  solution. 

In  regions  1-3  the  distance  to  the  magnetic  axis  is  much  smaller 
than  the  distance  to  the  symmetry  axis,  allowing  us  to  use  a  large 
aspect  ratio  approximation  there.  Using  Eq.  (4)  to  specify  the  poloidal 
current  distribution,  the  large  aspect  ratio  limit  of  the  Grad-Shaf ranov 
equation  gives 


3v<v$>+*2i3v<i§>-0  •  (A1> 

2 

where  V  =  r  ,  and  r  is  the  distance  to  the  magnetic  axis.  In  obtaining 
this  equation  we  have  specialized  to  circular  flux  surfaces  near  the 
magnetic  axis. 

In  region  l  we  expand 

\l>  ■  +  iJijV  +  +2^  +  •••*  •  (A2) 

Taylor  expanding  the  given  q(ip)  ■  qQ  +  (3q/3^)  +  ...,we 

substitute  it  along  with  Eq.  (A2)  into  Eq.  (Al)  to  determine 
i>2>  ♦3*  etc.  *n  terms  of  i|>o  and  1^.  We  get  for  <(>2» 

*2  *  "*1(1  +  qo  * 
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Its  contribution  to  i|>  is  small  everywhere  in  region  1  if  ij^a  «  <|>j» 

2 

where  6a  =  V(*  -6ib  ).  This  condition  will  be  satisfied  self- 
To  To 

consistently  by  our  solution.  The  constants  and  are  to  be 
determined  by  matching. 

Region  2  is  assumed  to  be  sufficiently  small  that  \|»  is 
approximately  constant  across  this  region.  In  addition,  we  integrate 
Eq.  (Al)  across  this  region  to  get  a  jump  condition  on 

♦  1  [6a2  +  R2(qQ-  6*fl  ^p)2]  -  +;  (6a2  +  R26q2/2)  ,  (A3) 

where  6q  is  the  value  of  6q  at  r**o. 
o 

In  region  (3)  we  expand  out  Eq.  (Al )  and  neglect  the  term 
containing  dq/dV.  Again,  this  approximation  is  justified  self- 
consistently,  taking  d6q/di|»  ~  6q/i|>.  The  solution  to  the  resulting 
equation  is 

♦  -  C.in  (r2  +  6q2  R2)  +  C0  ,  (A4) 

where  Cj  and  C2  are  constants  to  be  determined  by  matching. 

In  region  4  we  can  neglect  6q  entirely,  getting  the  standard 
elliptic  Integral  solution  for  the  field  of  a  ring  with  current  I  and 
major  radius  RQ, 

*  -  (4IR  /ck2)(R2-Hl2+2RR  sin0)”l/2[(2-k2)K(k2)-2E(k2)]  (A5a) 

000 

k2  -  4ftR  (ft2  +  R2  +2RR  sine)"1  ,  (A5b) 

000* 


where  R,0  are  spherical  coordinates  relative  to  the  center  of  the  ring, 
and  E  and  K  are  complete  elliptic  integrals  (we  follow  here  the  notation 

on 

of  Abramowitz  and  Stegun  ). 

Matching  all  our  solutions  together,  we  get  the  transcendental 
equation. 

y  »  Hln[6q2(ay+1)]  , 

where 

y  =  a|»o/I  , 

and 

a  =  •  (A6c) 

This  gives 

I  «  a|»  /in(l/5q2)  •  (A7) 

o  o 

Expressing  Air  j/c  *  XB,  we  have  X  =  2/(RqQ)  in  region  l,  and 
X  *  2MqQ(r2  +  6q2R2)  ln  region  3. 

Now  consider  the  equilibrium  specified  by  Eqs.  (13)  and  (14).  We 
again  use  the  same  approach,  hut  now  we  need  to  solve  separately  In  only 
two  different  regions.  Near  the  magnetic  axis  we  use  the  same 
approximation  used  in  region  3  of  the  previous  equilibrium.  Far  from 
the  magnetic  axis  we  use  the  elliptic  integral  solution  used  in  region  4 
previously.  Matching  these  solutions,  we  again  get  Eq.  (A7). 


(A6a) 

(A6b) 
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Fig.  7. 


Flux  surfaces  during  axisymmetric  coalescence 

a)  Initial  approach  of  the  two  spheromaks 

b)  Reconnection  of  the  magnetic  islands 

c)  Final  equilibrium  after  coalescence. 

Axisymmetric  coalescence  of  spheromaks  with  different  values 
of  ij>  on  the  magnetic  axis,  and  On  the  left  we  show 

the  q  profiles  before  coalescence,  and  on  the  right  we  have 
q(\p)  after  coalescence. 

Restoring  the  q  profile  by  purely  axisymmetric  coalescence. 

a)  q(iji)  of  the  fusion  spheromak  at  t-0  and  after  a  small 
time  has  elapsed 

b)  Refueling  profiles 

c)  q(ij>)  after  axisymmetric  reconnection. 

The  q  profiles  before  and  after  axisrnmetric  coalescence  for 
the  refueling  scheme  which  makes  use  of  an  n*l,  m*l  tearing 
mode. 

The  n»l,  m*l  helical  flux  function  before  (solid  line)  and 
after  (dashed  line)  the  growth  of  the  n*l,  m*l  tearing  mode. 
The  q  profiles  before  and  after  axisymmetric  coalescence  for 
a  refueling  scheme  which  uses  an  n*l,  m»l  double-tearing 
mode. 

The  n*l,  m-1  helical  flux  function  before  (solid  line)  and 
after  (dashed  line)  the  growth  of  the  n«*l,  m*l  double-tearing 
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